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Abstract 



^ ' Fission fragments from 1 A GeV ^^^U projectiles irradiating a hydrogen target were 

investigated by using the fragment separator FRS for magnetic selection of reaction 
CJ \ products including ray-tracing and AE-ToF techniques. The momentum spectra of 



733 identified fragments were analysed to provide isotopic production cross sections, 
fission-fragment velocities and recoil momenta of the fissioning parent nuclei. Be- 
sides their general relevance, these quantities are also demanded for applications. 
^ ' Calculations and simulations with codes commonly used and recently developed or 

. improved are compared to the data. 
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1 Introduction 



Fission of heavy nuclei excited at high coUision energies has been investigated 
for the last decades. Produced by fragmentation or proton collision, the ex- 
cited nuclei release their energy by emitting neutrons and charged particles, 
and eventually by fission. At low excitation energy, close to the fission barrier 
the probabihty of fission for nuclei around U stays constant or decreases at 
increasing excitation energy (first chance fission). Increasing further the exci- 
tation energy, two opposite tcndancics arc observed. High excitation energy 
and emission of charged particles reduce the probability of fission. On the 
other hand, neutron emission produces more and more fissile nuclei, enhanc- 
ing fission. For the highest excitation energies, lighter elements which are less 
fissile, are created and finally fission vanishes. The competition between the 
various decay channels provides keys to understand the structure of hot nuclei, 
level densities, dissipation and disappearance of shell effects with excitation 
energy. 

One of the first experimental approaches to the question of competition of 
the different decay-channels was opened with radiochemical and off-line mass 
separator techniques. G. Fricdlander et al. [1] studied isotopic production of 
the alkaline elements rubidium and cesium in p -|- U and p + W reactions at 
various energies. This pioneering work already characterises the evolution of 
the symmetric breaking of the excited fragmentation residues as a function of 
the incident energy. A decade later, similar measurements on alkalines were 
undertaken using on-line mass separators [2,3] which we will use as a basis for 
comparison in the following. Meanwhile many results evaluated in reference 
[4] were obtained to complete these first studies. 

The energy dependence of the reactions already studied in reference [1] must 
be known and controlled for all applications. Until now this aspect has been 
explored by radiochemistry and 7-spcctroscopy. The systems p -|- Pb [5] and 
p -|- Bi [6] were extensively investigated by using off-line 7-spectrometry in a 
wide range of energies. 

Why do we start new studies on this question ? The new experimental equip- 
ment operating at GSI offers an efficient way to separate isotopes in inverse 
kinematics at relativistic energies. A 1 A GeV U-beam is delivered at the 
heavy-ion synchrotron (SIS), and in-flight separation techniques become ap- 
plicable to all the energetic U residues. The inverse kinematics at relativistic 
energies provides fully ionised forward-focussed fragments which are separated 
by using the FRagment Separator (FRS) and identified by energy loss (AE) 
and time-of-fiight (ToF) measurements. Each of the isotopic cross-sections can 
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be determined, and the kinematics of every residue can also be reconstructed. 
Based only on physical properties of the ions, the method does not depend on 
chemistry. As beta-decay half-lives are much longer than the 0.3 /is separation 
time, primary residues are observed. 

One of the important findings of the combination of the relativistic U beam 
at SIS with the FRS separation techniques was the identification and yield 
measurement of 117 new neutron-rich nuclei down to the sub-nanobarn level, 
among which the doubly- magic nucleus ^^Ni [7,8]. The efficiency, selectivity 
and high sensitivity of the method was clearly demonstrated. It triggered new 
efforts to produce beams of exotic nuclei by in-flight techniques. 

Apart from the scientific interest of proton induced spallation reactions, appli- 
cations were soon looked upon. For the purpose of Accelerator Driven Systems 
(ADS), for energy production and transmutation of nuclear waste [9], a great 
number of data sets were missing. Spallation targets for neutron-sources re- 
quire extensive data in the same domain. Patents were approved in Germany 
and France [10]. Thus a campaign of systematic measurements was under- 
taken to obtain isotopic production cross-sections for a number of systems 
with an accuracy of about 10% and an improved understanding of reaction 
mechanisms in order to interpolate to other systems and energies below 1 A 
GeV. 

For this purpose, at 800 A MeV Au + p [11,12,13], Pb + p at 1 A GeV [14,15] 
and at 0.5 A GeV, Pb + d at 1 A GeV [16] and spallation from 1 A GeV U 
-I- p [17] and 1 A GeV U -I- d [18] were already studied and pubhshed. The 
missing part is still under analysis. Spallation of Fe on p and d was measured 
between 0.3 and 1.5 A GeV. Very recently, the systems Xe + p and Xe + 
d were measured between 0.2 and 1 A GeV to conclude this experimental 
program. 

In this article we report on the (U + p)-fission products measured at 1 A GeV. 
The system plays the role of a prototype for high-energy nucleon-induced fis- 
sion of an actinide. At high energies, data on fission of ^^^U induced by protons 
can be used as a benchmark for other actinide nuclei present in incineration 
situations. 

The results on fission-residues in spallation of 1 A GeV ^^^U on protons to- 
gether with our results on evaporation-residues [17] bring a comprehensive 
overview on this system. Earlier results on U-projectile fission were obtained 
on Pb [19,20,21,22] and on Be [23] targets. Altogether they should lead to 
a coherent reconstruction of the intermediate systems after the primary cas- 
cade and nucleon evaporation preceeding fission. The reproduction of the data 
with standard codes, the intra-nuclear cascade-evaporation model LAHET 
[24] is far from being satisfactory. Improvements of codes have been published 
[25,26,27,28], and are applied here. 
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Finally, the new data on isotopic yield measurements for ^^^U and ^°^Pb to- 
gether with present and coming data on their energy dependence should pro- 
vide a reference in order to control nuclide production in the forthcoming 
applications of proton-induced spallation. 



2 Experiment 

2.1 Experimental set-up 

Uranium-ions delivered by the Penning source are accelerated in the UNILAC 
and post-accelerated in the SIS of GSI-Darmstadt. The beam energy of the 
238^76+ jQj^g target is 1 A GeV and their velocity — 0.876. The average 
time structure is a pulse of 7 s every 13 s. The beam intensity delivered is 10^ 
to 10^ ions/pulse. It is controlled and monitored prior to the target by means 
of a secondary electron chamber, called SEETRAM [29] . 

A plastic target could not be chosen, since the production yields from C prevail 
upon the H yields in the mass region intermediate between fragmentation and 
fission. The uncertainties in the sharing between C and H in a plastic-target 
would not be consistent with our accuracy requirements. 

The liquid hydrogen target used in the present experiment consist of a (1.25 ± 
0.03) cm thick circular cell closed by two thin parallel titanium windows. The 
cell is located at the bottom of an helium cryostat. Surrounding the cell, thin 
aluminised mylar reflectors and external titanium windows ensure thermal and 
vacuum insulation [30]. The hydrogen thickness of (87.2 ±2.2) mg/cm^ [11] is 
chosen to obtain a 10% reaction probability. The Ti target walls have a total 
thickness of (36.3) mg/cm^. Their contribution to the counting rates remains 
below 10%. It has been measured systematically using a "dummy" target of 
the same composition as the empty target cell. 



2.2 Separation and identification in A and Z 

At energies close to 1 A GeV, fission fragments are observed up to Z = 72 
and they are nearly totally stripped. In the worst case at Z = 72, less than 
15% of the ions still carry electrons [31]. They are momentum analysed with 
the fragment separator FRS, a double-stage achromatic spectrometer with an 
intermediate dispersive image plane S2 and a final image plane S4 [32]. A 
schematic view of the FRS and detectors is shown in Fig. 1. 
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Fig. 1. Schematic view of the fragment separator with the detector equipment. 

In S4 the separated fragments are identified in Z by their energy losses mea- 
sured in a four-stage ionization chamber, MUSIC [33]. The mass number A 
of the analysed fragments is determined from a measurement of velocity and 
magnetic rigidity between S2 and S4. The measurement is performed by two 
position sensitive plastic scintillators [34]. The times and positions in S2 and 
S4 have to be known accurately in order to correlate the time of flight with the 
precise length of the flight-path. Fission fragments fill the full phase space de- 
fined by the acceptance of the FRS. The presence of the 5 mm thick scintillator 
in S2 slightly alters the achromaticity of the FRS. Nevertheless 36 elements 
from Z = 28 to 64 are simultaneously transmitted to S4. Signals from the 
S4 scintillator trigger the data-acquisition system. Mass and charge numbers 
are calibrated by comparison to the parameters of the primary uranium beam 
and by using the known structure of fission yields in the asymmetric fission 
domain. 

2.3 Constitution of the momentum spectra 

The momentum acceptance of the FRS is 3 %. Applying the basic relation 



with Z and A numbers determined, a 3 % fraction of the (3j spectrum is 
obtained for each measurement. The total spectrum of transmitted fission 
fragments covers a range of 7 % to 20% of the reduced momentum of the 
projectile, /9o7o- The measurements consist of a 2 % step-by-step scanning of 
the FRS magnetic fields tuned on cadmium, (Z = 48). After normalisation 
to the incoming dose given by the SEETRAM, and dead-time corrections, 
complete (3j distributions are reconstructed by assembling about 10 files of 
overlapping Bp settings to produce a (3^ spectrum for each of the isotopes; 
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where iij (/97) is the i^^ bin of the spectrum, Tj is the dead time fraction 
and Nj the number of incident ions for this measurement. All in all 4-10^ events 
were accumulated. 

Fragments arc spread over 20 cm in the image plane at S2 covered by the first 
position-sensitive scintillator. For an identified isotope the distribution in S2 
depends only upon /37, and the accuracy on S2 position governs the accuracy 
of the magnetic rigidity. In the present case, a precision in position of 2.5 mm 
(FWHM) leads to a resolving power of (3^ / Ap^ — 2600, given the dispersion 
in momentum of the first stage of the FRS by percent, D = 6.5 cm/%. In 
this experiment the response of the scintillator was slightly non-linear, and a 
routine was systematically used in the analysis software to correct for this, 
[11]. 

The FRS magnetic fields tuned on Z = 48 guide this element to the center of 
S4 and display on both sides the neighbouring elements, gradually apart from 
the center of S4. The 36 elements are transmitted together to the final focus 
S4, and momentum distributions of the produced isotopes arc simultaneously 
measured in one series of Bp settings of the FRS. The (5j spectra can be con- 
verted to the projectile frame and are presented in cm/ns as velocity spectra 
in the system of the projectile. This coordinate system is appropriate for com- 
parison with previous results concerning fission from target nuclei. Moreover, 
at high excitation energy, fission fragments are emitted isotropically in the 
fissioning system [19,35]. This allows for analytical formulations to be used in 
the analysis of the distributions. 



2.4 Description of the velocity spectra 

The reconstructed velocity spectra are the corner stones of our analysis. An 
example of the velocity spectrum is shown in Fig. 2b. Two narrow asym- 
metric peaks are observed roughly equidistant from velocity zero, i. c. the 
projectile velocity. Velocities are defined in reference to the beam direction 
i. e. forwardly and backwardly emitted fragments have positive and negative 
velocities, respectively. The velocities of the fission fragments arise from their 
Coulomb repulsion at scission. Velocity vectors of a specific isotope populate 
a thin spherical shell in the fissioning system, which may be slightly shifted 
compared to the projectile system by the primary reaction recoil. The sphere 
is pictured in Fig. 2a by a circle, the cut of the sphere by a plane which con- 
tains the beam. Only forward and backward cups of the sphere, defined by 
the angular acceptance of the FRS a, are transmitted, and the longitudinal 
projections of their velocity distributions are shaping the two peaks (Fig. 2b). 
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Fig. 2. a) Schematic view of the experimental parameters shaping the measured 
velocity spectrum in the frame of the fissioning system. V/ is the fission-fragment 
velocity, a is the angular acceptance of the FRS, and cJq, its variance, tp^'^ are 
the corresponding emission angles of fission fragments in the fissioning system in 
forward and backward directions, respectively, b) Velocity spectrum of ^^^Te in the 
frame of the fissioning system. The velocity V = refers to the projectile frame. 
Vapp is the apparent fission velocity defined in the text. 

For selected elements the velocity spectra of all isotopes are presented on 
scatter-plots in Fig. 3. Rates are given in terms of velocities and of neutron 
number N. They all show strongly populated areas along two parallel lines 
associated to backward and forward emitted fission fragments. The lines, al- 
most symmetric in respect to the axis (V = 0, i.e. projectile frame), show 
a fission fragment velocity decreasing with increasing atomic number Z of 
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the fragment, a trend due to momentum conservation in fission. The parallel 
lines contrast with U + Pb scatter-plots [22] where 1) neutron-rich nuclei are 
strongly enhanced for zirconium and tellurium due to electromagnetic fission 
and 2) fragmentation residues sit in between the forward and backward regions 
for neutron-deficient isotopes. In spite of the absence of evaporation residues, 
at decreasing masses of the lighter elements, events with small velocities are 
found. Besides the small contribution of fragmentation residues produced in 
target windows, this is a first indication for the contribution of secondary re- 
actions which will be discussed in the following section. 




Fig. 3. Cluster plots of number of events measured on target in terms of velocity 
versus neutron number. The fission velocity coordinate is given in the center of 
mass of the projectile. The intensity scale is logarithmic with a factor of 2 between 
adjacent colours. 

Data taken on the "dummy" target show the same pattern as in Fig. 3 with 
a contribution of fragmentation. Yields are defined by the integrated velocity 
spectra. Isotopic distributions of yields in the "dummy" and in the target are 
compared in Fig. 4. As intended, the yields in the "dummy" are only a few % 
of the yields in the target, except in the neutron-deficient region where they 
reach larger values, but less than 20% . Asymmetric fission, leading to neutron 
rich isotopes, is also clearly enhanced for titanium relatively to hydrogen. The 
hydrogen contribution is obtained by subtracting the yields on the "dummy" 
from the target yields. Uncertainties are estimated from the fluctuations of 
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the intensity ratio of the two peaks, as will be discussed later. The error bars 
on yields measured for the "dummy" are negligible when reported to yields 
on hydrogen. 




N N N 



Fig. 4. Comparison of transmitted yields [arb. units] of isotopes on hydrogen 
(squares) and "dummy" (circles) targets as a function of their neutron number 
for selected elements. 

Evaluating velocities, the /?7 distributions on target and "dummy" are found 
to be very similar for a given isotope, and the yield of the "dummy" is on a 
few percent level. Therefore, the contribution of the titanium windows could 
be ignored in the evaluation of fission-fragment velocities in U + p. 

2.5 Fission velocity and transmission 

The geometrical aperture of the spectrometer seen from the target defines 
a solid angle fl simply related to the angular acceptance a assuming axial 
symmetry, a is the limiting value for the emission angle of transmitted fission 
fragments in the laboratory system. corresponds to angles ip^ and ip^ in the 
fissioning system with F and B set for forward and backward, respectively, 
see Fig. 2a. The angular transmission T^ is defined as the ratio of the yield 
measured for a given fission fragment to the total production yield for this 
fragment. 

The determination of Tq is a key to obtain the cross-sections. It depends 
upon fission fragment velocities and it must be evaluated for each fragment. 
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Tq is a function of the three variables a, Po representing the beam and (3f 
the fission process with V/ = /3/.c, the fission velocity, i.e. the radius of the 
sphere shown in Fig. 2a. As shown in Appendix 1, for a bundle of particles 
with an angular divergence a starting from a point-like source and entering 
the FRS on its principal trajectory, the transmission Tq can be calculated 
rigourously, without any approximation, when setting the variances to zero. 
In the fissioning system, Tq is obtained by integrating the surface elements of 
the spherical velocity shell within the two solid angles, and we obtain [22]: 



_ ^ o ^ p cos f/?-'^ + cos </7-^ 

Tn = Tn^ + Tn^ = 1 ^— ^ (3) 



In Fig. 2b the velocity distribution of transmitted ^^^Te is given as an exam- 
ple. Three quantities are evaluated from these distributions: 

1) The spectrum dn{(3'~f) / d{f3j) integrated over the whole range of P^, eval- 
uated for each of the isotopes gives the measured yields. 

2) A shift of the mean velocity of forward and backward emitted fission frag- 
ments, compared to the projectile velocity (V = 0) gives the velocity of the 
fissioning parent nucleus, the recoil velocity. 

3) The centroid of each of the two peaks defines a mean fission-fragment ve- 
locity, the apparent velocity = Vapp/c different for forward and backward 
emitted fragments. We define the mean value of the two apparent velocities. 



Papp = ' ' I ' ' (4) 

The two peaks are shaped by the variances of the three variables a, /3o and 
(3f, as will be shown later by a Monte-Carlo simulation. With all variances set 
to zero -Fig. 5a- the peaks become rectangles, as the phase-space density of a 
spherical homogeneously filled shell per interval d/3\\ with the longitudinal 
velocity, is constant. The projection of the spherical shell on the beam axis is a 
rectangle. In this case the apparent velocity of forward and backward emitted 
fission fragments (see Fig. 2) is given by: 



PB a Pf - PfCOSifi^'^ 1 + cosy?^^'^ 
Papp ^Pf 5 = Pf 5 



and the mean value follows: 



&pp=^(i+ ^°^'^'r°^^° ) (5) 
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With eq. (3) defining the transmission, we obtain the fission velocity depending 
on (3app and Tq: 



2-Tn ^ ' 



from which the limiting cases are deduced /?/ = (3app for Tq = 0, and (3f — 
2f5app for Tn = 1. 

In Appendix 1 we show that with ip-^ — {(fi^ + </?^)/2 the fission fragment 
velocity /3/ and the transmission Tq are given in very good approximation by: 



/5/ = /3app[l + (tan(^s/2)'] (7) 



- l + [cot^./2]2 



with 



tan(^s/2) = ^ (9) 

■^Papp 

Piini is defined by a and the velocity of the fissioning parent nucleus /3 approx- 
imated by Po 



Aim = aV7^ - 1 ~ Q;/?o7o (10) 

Eq.(7) and (8) demonstrate the central importance of the ratio Pum/'^Papp- 
Summarising the procedure: 

• Prom the distribution, the mean apparent velocity Papp — {Papp + f^app) /2 
is evaluated for each isotope. 

• Tq and Pf are simple functions of the mean opening angle ips = {<P^ + 

• Aim is related to the acceptance angle a of the spectrometer and to the 
projectile velocity (5q 
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2. 6 Simulation of the angular cut of the FRS 

In order to include the variances on the FRS angular acceptance and the 
kinematics of fissioning nuclei, we have simulated the physical cut in phase 
space due to the angular acceptance of the FRS. The purpose is to relate V/ 
and the transmission to the measured Yapp- 

The Monte-Carlo calculation is based on the following assumptions: 

1) In the system of the fissioning parent nucleus, approximated by the system 
of the projectile nucleus, fission fragments are isotropically emitted with a 
given velocity V/. 

2) The spherical shell of velocity vectors is converted to the laboratory system 
moving relativistically with a velocity (3. Fragments emitted within the FRS 
acceptance angle a are transmitted. At this stage we obtain a transmission 
value and the momentum spectrum in the laboratory system. 

3) This spectrum is reconverted into the projectile system. 

4) Variances of velocities and of the acceptance angle are further parameters 
which are included. 

The simulation enlightens the effects of the experimental constraints on the 
measurement. On Fig. 5a the forward and backward rectangles correspond 
to the share of transmitted fission- fragments with ctq = 0. A filled rectangle 
would be obtained for j3f smaller than af3'-f and T^ would then reach 100 %. 
The second curve in the same frame results from the introduction of a variance 
of the acceptance angle cr^ = 2.5 mr derived from the known geometry of the 
FRS [36,37]. It shows that the filling of the spectra at intermediate velocities 
comes mainly from the large variance of the FRS acceptance angle. 

Fig. 2b versus Fig. 2a shows that the external slopes depend only on the fiuctu- 
ations of the reduced momenta of the fission fragments and/or of the fissioning 
nuclei. Several sources contribute to these fiuctuations - fission fragment ve- 
locities, target location straggling and recoil momenta of the fissioning nuclei 
-, and lead to a broadening of the peaks as illustrated on Fig. 5b. The relative 
contributions of these sources are quantified in the discussion, sec sect. 4.8. 

With a/s-y = 7.10"^ the external slopes are almost vertical and the value of V/ 
= 1.2 cm/ns introduced in the simulation can be extracted from the distance 
between the external sides at half maximum. On the contrary, if a^^ becomes 
larger than the width of the rectangles of Fig. 5a, this half-maximum rule is no 
longer valid, as seen for a^j — 14.10"^. A value of ajjj — 7.10"^ is selected as 
an "empirical broadening" to optimize the simulation. The result is illustrated 
in both examples of Fig. 6a and 6b. 

The relation between Yapp ^-nd V/ either from eq. 6 or from the simulation 
are very close as shown in Fig 7(a) and the trends at Tq = or 1 are correct 
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Fig. 5. Impact of variances a) ctq, of a, the acceptance angle of the FRS, and b) 
ap^ of (3^, the momentum of the fissioning nucleus. Fission-velocities are different 
in a) and b) in order to better illustrate the contributions of both variances. 




-2-1012 -2-1012 
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Fig. 6. Examples of simulation with ctq, =2.5 mr ap^ = 7.10 ^ compared with the 
measured spectra for ^^^Te and ^^''Sm. 

in both cases. Simulated and calculated values of Tq are the same for fission 
velocities larger than 0.8 cm/ns ( Fig. 7b). For smaller velocities the calcula- 
tions, which do not include the large variance of the acceptance angle, cannot 
describe the transmission. Approaching Tq = 1 the simulation accounting for 
the variances provides better results. It is seen on the velocity spectra of iso- 
topes of neodymium to gadolinium where a few isotopes (the heaviest isotopes 
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Fig. 7. Fission fragment velocity as a function of apparent velocity a) and trans- 
mission b) as a function of fission fragment velocity given by the simulation (full 
points) and calculations (empty points) using eqs. (7) and (8). Lines are fits to the 
full points. 

showing the largest fission velocities) still show a residual valley between for- 
ward and backward peaks. In the following transmissions and fission fragment 
velocities are obtained using the values of N app determined for each isotope 
and the fitted curves of Fig. 7. 



3 Experimental results 

3. 1 Fission fragment velocities 

The velocity spectrum of each isotope has been analysed to extract a value 
of the apparent velocity. The apparent velocities have been converted into 
fission fragment velocities. For a given element, the fission fragment velocity 
decreases smoothly from a maximum value for the heaviest isotopes down to 
values smaller by 10% for the lighter isotopes, Fig. 8. Further down to the most 
neutron-deficient isotopes, the velocities fall abruptly. In case of U+p, this fall 
cannot be related to the occurence of evaporation residues since this process 
does not populate the neutron- deficient region of fission fragments. The small 
contribution from fragmentation yields in the Ti-windows (Fig. 4) cannot ex- 
plain the observed fall in velocity of 20%. This fall is due to secondary reactions 
in which heavier primary fission fragments loose nucleons in a secondary frag- 
mentation reaction in the target. Atomic numbers and masses are reduced, 
whereas the velocity spectra are only blurred. These contributions - due to 
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Fig. 8. Fission velocities for selected elements. Cross symbols correspond to regions 
where secondary-reaction products are mixed with direct fission fragments. 

larger Z with smaller fission-fragment velocities - are superimposed with the 
spectrum of the primary isotope. This contamination reduces the apparent 
velocity of the isotope. It becomes significant in the neutron-deficient region 
where productions by direct fission is low and secondary-reaction products are 
abundant. However, the small velocities of the lightest isotopes of neodymium 
cannot be explained by this argument, as secondary reactions of fission frag- 
ments contribute little. Here secondary reactions of spallation products could 
pollute the yields of neutron-deficient fission fragments. 

The velocity of an element is taken as the mean value of the velocities for 
the four isotopes most abundantly populated, in order to select the dominant 
symmetric fission regime. This velocity is plotted as a function of Z in Fig. 
9a. The uncertainties are less than 4%. Three fines are indicated showing 
the velocity dependence expected for three values of the nuclear charge of the 
fissioning parent nucleus Zq = 92, 86 and 80 [38,39]. The scission configuration 
of the fissioning nuclei is approximated by the outline of two deformed nuclei 
{ Pi = (^2 = 0.65) at a given distance. Using momentum conservation, the 
fragment velocities are calculated for the total kinetic energies released from 
the scission configuration. 

None of the three lines coincides exactly with the data. Fissioning nuclei 
with Zq between 84 and 90 can contribute to the production of fission frag- 
ments. The general trend towards smaller fission fragment velocities going 
from neutron-rich isotopes to neutron-deficient isotopes is observed for all el- 
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Fig. 9. Fission fragment velocities as a function of Z. a) The lines are obtained using 
the calculation of velocity for fissioning parent nuclei of Zq = 80 (full line), Zq = 86 
(dashed line) and Zq = 92 (dotted line) with the parametrisation of [38] [39] and 
a deformation of /3 = 0.65 for both fragments, b) Measured variances of the fission 
velocities in cm/ns. c) Measured displacement of the mean velocity of forward and 
backward emitted fragments. The lines indicate the effective broadening b) and the 
recoil velocity c) as measured for the mean parent-fissioning nucleus [17]- see text 
sect. 4.8 -. 



ements (Fig. 8) corroborating that heavy (Z = 90) and light (Z = 84) parent 
nuclei do contribute. The first ones contribute to the more neutron-rich and 
the last to less neutron-rich isotopes. 
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The slope of the external edges of the velocity distributions introduced empir- 
ically in the simulation (Fig. 5b) as a variance up-y are now evaluated system- 
atically as a function of Z and plotted in fig. 9b). Unfolding these measured 
values and the underlying rectangle of the transmitted momenta, fig. 5a, the 
effective broadening ay^ is obtained, and given by the line in fig. 9b. Both 
variances are related by cry^ = ^ • (7/3-y • c. Within the limit of our accuracy 
this width does not depend on the atomic number. Its mean value is — 
0.13 cm/ns, which corresponds to ag^ ^ 8.10"^ in good agreement with the 
value of 7.10"^ empirically deduced in sect. 2.6, from the simulations of the 
measured spectra (Fig. 5b). 

The recoil of the fissioning nucleus is determined from the mean value of 
the forward and backward fission velocities. The values plotted in Fig. 9c are 
averaged over the four most abundant isotopes of an element. This recoil is 
small and it decreases from -0.03 at Z = 62 to -0.13 cm/ns for Z = 30. We 
obtain a mean value of -0.08 cm/ns. 

3.2 From yields to cross sections 

Isotopic cross sections of (U + p)-fission are obtained by dividing the produc- 
tion yields on hydrogen by the proper transmission, as presented in sections 
(2.5) and (2.6). Our final values arc compared to those used in the analysis 
of the parallel experiment U -|- d at 1 A GeV evaluated by J. Pereira apply- 
ing a different method [36]. A good agreement is found between both sets of 
transmission- values evaluated in terms of the fission-fragment velocities. 

The correction for secondary reactions occuring in the target has to be applied. 
All fission fragments, neutron-rich more than proton-rich undergo secondary 
spallation and the yield is depleted. For neutron-deficient isotopes, the fis- 
sion yields are increased by secondary evaporation residues. The correction 
for secondary reactions enhances the share of neutron rich isotopes and re- 
duces neutron-deficient ones. This is the reason for the "hook" shape of this 
correction already presented in a previous article [15] and calculated more pre- 
cisely in a forthcoming paper [41]. When they exceed 50 %, cross sections for 
corresponding isotopes are no longer considered. Our experimental method to 
determine cross sections of primary residues reaches its limit. 

Moreover, small secondary effects need to be corrected for: For fragments 
not centered in S4, i. e. apart from the optical axis, the FRS acceptance 
angle and consequently the transmission are slightly reduced. The FRS angular 
acceptance falls from 14.8 mr on the optical axis down to 13.3 mr at 8 cm apart 
from the axis in S4 [36] . A smooth parabolic function was used to account for 
this relative loss in transmission. The correction given in table 1, reaches 14% 
for elements like nickel and samarium. 
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Table 1 



Correction factors appHed to the measured yields to account for losses due to 1) 
the shift in S4 position compared to the central position, 2) the secondary reactions 
in the SC2-scintillator at S2, and 3) the atomic charge exchange in target and 
S2-scintillator for Z > 52. 

Fragments reacting in the scintillator SC2 at S2 escape the magnetic filter of 
the second half of the FRS. The related loss increases hnearly from 11 % for 
nickel up to 17 % for gadolinium. It is calculated using the total reaction cross 
section formulation of Karol [42] . The results are given in table 1 for the range 
of elements investigated. 

A loss by ionic charge exchange occurs, either in the target or in the plastic 
scintillator S2. Eventhough the probability is small for fission fragments at 
relativistic energies, this factor, negligible up to tellurium, reaches 9 % for 
gadolinium. Values, calculated using the code GLOBAL [31] are to be consid- 
ered for heavy elements as shown in table 1. 

3.3 Discussion of uncertainties 
a) Systematic uncertainties 

• The uncertainty on the SEETRAM calibration is evaluated to 7 % [29] 

• The uncertainty on the transmission is evaluated to 5 %. It is mainly due 
to the 4 % uncertainty on the fission-fragment velocity. 

• The uncertainty on the target thickness was established in a dedicated ex- 
periment to 2.5 % [11] for a beam hitting the target in its center. However, 
an uncertainty of 6 mm in the vertical position of the cryostat supporting 
the H2 cell leads to an increased uncertainty of 4 % on the target thickness 
[17]. 
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• The loss factors given in table 1 are determined with an overall uncertainty 
of 3 %. 

All contributions quadratically added up lead to a total systematic uncertainty 
of 10 %. This systematic uncertainty is not included in the error bars of the 
cross sections presented in tables of Appendix 2 and in figures 10-12. 

b) Statistical uncertainties 

• The single-event uncertainty on the angular transmission is large and mainly 

determined by the large variance of a. a a is evaluated to 17 % [36]. In com- 
parison the uncertainty on (iapp is negligible. Taking 50 observed events 
corresponding to the smallest cross section (20/ib) as the worst case, the 
statistical uncertainty for the transmission reaches 4 %. Statistical uncer- 
tainties are smaller than the systematic uncertainty. 

• Uncertainties for the yields have been determined from the ratio of the 
rates integrated for forward and backward emitted fragments. The difference 
between the measured ratio and the kinematically expected ratio, is taken 
as twice the relative uncertainty of the yields. The absolute uncertainties 
are given in Fig. 4, Fig. 10-12 and in Appendix 2. For the most abundant 
isotopes the relative uncertainty is approximately 2 %, that is much less than 
the systematic uncertainty. For the less abundant isotopes, the statistical 
uncertainties are dominant and overcome the total systematic uncertainty. 

Isotopic cross sections are given in Appendix 2 and their distributions for each 
of the elements appear in Fig. 10-12. They show a rather regular bell shape 
with some enhancement in the region of neutron-rich isotopes for elements 
known as asymmetric fission products. At the maximum of the elemental dis- 
tributions, cross sections vary from 0.5 mb to 20 mb. Values down to 20 //b 
are measured, mostly in the neutron-rich region not polluted by secondary re- 
actions. On the neutron-deficient side, high rates of secondary reactions limit 
cross-section measurements to values larger than a few mb. 
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Fig. 10. Isotopic cross sections for fission-fragments in the reaction U 1 A GeV 
+ p for elements between 28Ni and 39 Y. 
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Fig. 11. Isotopic cross sections for fission-fragments in the reaction U 1 A GeV 
+ p for elements between 40 Zr and 51 Sb. 
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Fig. 12. Isotopic cross sections for fission-fragments in the reaction U 1 A GeV 
+ p for elements between 52Te and esEu. 
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3-4 Comparison with previous isotopic distributions 



Measurements of primary isotopic cross sections had been performed for alka- 
lines only. Using fast on-line separator techniques combined with the chemical 
selectivity of surface ionisation ion-sources [2,3] relative cross sections were 
obtained covering a wide range of isotopes. The normalisation of the yields 
to cross sections was based on the results of G. Fricdlander [1]. Production 
yields on U-fission were measured at a proton energy of 24 GeV for isotopes 
of the complementary elements rubidium and cesium [43]. The distributions 
are compared to our results in Fig. 13. At 24 GeV, besides fission, evaporation 
residues populate all elements down to the region covered by fission. Moreover, 
the excitation energy range covers a wider domain, and evaporation of charged 
particles and neutrons leads to more neutron-deficient isotopes. The elemen- 
tal cross sections for both rubidium and cesium are about 75 % of the values 
presently measured at 1 GeV, in the same ratio as total fission cross sections 
[4]. In the region of the 8 most neutron-rich fragments, isotopic cross sections 
are the same as the ones presently measured. The probability to excite the U 
nucleus with a small excitation energy to produce neutron-rich fragments is 
independent of the proton energy, as indeed mentioned already in [1]. 
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Fig. 13. Comparison of isotopic distributions previously measured by J. Chaumont 
[43] -dotted-dashed line- and by B. Belyaev [3] -dashed line- with the present data 
-full points-. 

Later the production of the isotopes of cesium and rubidium was investigated 
for the system ^^^U -|- p at 1 GeV by using also on-fine mass separator tech- 
niques [3]. For rubidium our data are very consistent, but cross sections are 
found larger by a factor 1.8 for the six cesium-isotopes with N>85, as seen 
in Fig. 13. In the present experiment, isotopic yields are measured simulta- 
neously for the 36 elements populated in fission, insuring a coherent mass 
calibration for cesium as for all other elements including rubidium. Moreover 
the extraction of the asymmetric fission component of all elements achieved 
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in sec. 4-3 confirms the present mass cafibration. Finally the agreement of our 
isotopic distributions with those of the light elements analysed in the work 
of [44], also validates our results for the cesium isotopic distribution. In the 
on-line mass-separator measurement [3] the separation efiiciency might have 
overestimated for these isotopes. 



4 Features of fission 

4-1 Elemental cross sections 

The elemental cross sections, integrated over isotopes for each element are 
plotted in Fig. 14a. The distribution shows a Gaussian shape with a shoulder 
around Z = 52 revealing the contribution of asymmetric fission. The mean 
value of Z is Z = 45. It suggests that 90 Th is the mean element undergoing 
fission. The variance of the Z-distribution is 7 charge units. 

In Fig. 14b the ratio of the mean number of neutrons N/Z for each element is 
given. The largest mean number of neutrons per element and coherently the 
largest variances are found around Z = 55. A mean number of neutrons N — 
63 is deduced. 

The local variance a^, plotted on Fig. 14c, increases from 2.5 to almost 5. 
The presence of the asymmetric fission component is revealed by the maxima 
at the two peaks of asymmetric fission. 

4-2 Total cross sections 

The fission cross section summed over all isotopes is found to be (1.45 ± 
0.15) b. However, a few elements are not involved in the analysis, in the region 
above Z = 64 and below 28. Above Z = 64 the cross section contributing is 
evaluated to 7 mb [18]. Below Z = 28, the analysis of M. V. Ricciardi [44] gives 
a contribution of 72 mb. Thus the total fission cross section amounts to (1.53 ± 
0.15) b. This value compares well with previous results of direct measurements 
of fission cross sections (1.46 ± 0.07) b and (1.48 ± 0.06) b [35,45], respectively. 
A value of (1.52 ± 0.12) b was obtained in a recent measurement at GSI [46]. It 
is worth noticing that the present technique offers a high degree of selectivity 
and sensitivity to measure isotopic yields but a lower precision for total cross 
section measurements than simpler dedicated experiments. If the cross section 
for evaporation residues of 0.46 b, obtained by summing isotopic values found 
by J. Taieb in his analysis of the ^^^U 1 A GeV -|- p reaction [17] is added, 
a reaction cross section of (1.99 ± 0.17) b is obtained. It agrees with the 
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value of 1.96 b obtained by using the Glauber approach as described by P. 
J. Karol [42] with updated nuclear-density distributions [47], and with the 
INCL-calculation 1.94 b [28]. 



4-3 Symmetric and asymmetric fission 




Fig. 14. Integrated distributions of fission fragments from (^^^U + p) at 1 A GeV 
(full points), symmetric fission (empty points) and asymmetric fission (crosses), a) 
and d): Z-distribution, b) and e) mean neutron number divided by Z as a function 
of Z, and c) and f) variance in neutron number for a given Z. The short line in 
frame e) corresponds to P = -0.04 [48]. The thick line in frame b) and e) shows the 
position of the stable isotopes. 

The isotopic distributions of elements above Z = 50 can be decomposed into 
fragments arising from two components: high-excitation symmetric fission and 
low-excitation asymmetric fission. A decomposition into two Gaussian distri- 
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butions has been performed, as illustrated in Fig. 15. For the neutron- rich 
component the mean mass number Ap and the dispersion af^ are found to 
agree with the values known in asymmetric fission, see for example C. Don- 
zaud et al. [20]. For the light fission fragment group, the share of asymmetric 
fission does not exhibit so clearly and the known values of Ap and af^ [20] are 
used as further inputs. 




Fig. 15. Examples of the shares of symmetric (full lines) and asymmetric (dotted 
lines) fission processes. The sum, resulting from the least square fit is indicated. 



In spite of the small yields of the asymmetric component, cross sections for the 
groups of light and heavy fragments are obtained. Z-integrated cross sections 
are found the same for the two groups, as to be expected, and equal to (105 ± 
10) mb. For symmetric fission a cross section of (1.42 ± 0.15) b follows. The 
distribution of elemental cross sections becomes a symmetric curve slightly 
shifted (empty points in Fig. 14d compared to the total distribution on Fig. 
14a. The parameters N / Z and variances cr^ related to the two modes are 
clearly different. Their values are compatible with two mean primary fissioning 
nuclei : for the asymmetric mode ^^'*U and for the symmetric mode ^^^Th, with 
6 post-scission neutrons added. 

The N / Z values of the fragments are not constant , but a slow and regular 
increment of N / Z with the Z of the fragment is observed for the high ener- 
getic symmetric fission process, see Fig 14e. An electric polarisation = 
- 6Af°'' = A2 - Z2XA0/Z0 >0 for the fission fragment pairs A2, Ai is indi- 
cated. Beyond the regime of nuclear structure dominating asymmetric fission 
(E* > 40 MeV) a constant charge polarisability of P = -0.04 is expected [48]. 
P can be related to the slope of the / Z-dependence in the range of sym- 
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metric fission (Zo/2 = 45): 



-8^\Zl_ d{N/Z) 

The thin hne in Fig. 14e shows a slope calculated with P = -0.04. 

On Fig. 14f the variance cr^ in the region of symmetric fission increases with 
the Z-value of the elements. On the contrary for asymmetric fission, as known 
for a long time [49], a narrow distribution of isotopes is populated with 
a large neutron excess N /X . Beyond Z = 55 the values of N / Z and of 
decrease again and the valley of stability is crossed at Z = 58. For elements 
still heavier, fission fragments populate the neutron-deficient side of the valley. 



4-4 Mass distribution 
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Fig. 16. Mass-distributions of fission fragments for ^^^U (1 A GeV) -|- p. 

The mass distribution is obtained by summing isotopic cross sections for a 
value of A fixed, Fig. 16. Since isobaric yields do not depend upon /9-decay 
half-lives, mass distribution are accessible also by radiochemical methods and 
spectroscopic methods [5,6,50]. The distribution shows a wide Gaussian-hke 
shape slightly asymmetric: the slope is steeper on the side of light masses. 

The mean mass is A = 108.0 ± 0.3 and the variance of the distribution aA = 
(17.5 ± 0.5) a.m.u. Enhancements on both sides result from the asymmetric 
fission component. Since elements below Z — 28 are neglected, which con- 
tribute by about 5% to the total cross section [44], the value of A given is 



27 



reaction 


^sym 

(barn) 


(MeV) 


Z 

(a.ch.u.) 


N 
(a.ni.u.) 


(a.ch.u.) 


(a.m.u.) 


(a.ch.u.) 


U + p 
U + Pb 


1.42±0.2 
1.4±0.2 


76±3 
79±2 


45.0±0.1 
42.9±0.6 


63.0±0.15 
58.1±0.3 


7.0± 0.2 
6.9±0.7 


17.5±0.5 
17.2 ±1.7 


1.2±0.1 
1.3±0.3 



Table 2 

Parameters characterising fission of ^'^^U(l A Gev) + p in comparison to the sym- 
metric share of the fission process for ^ss^j _j_ pj-, Q.^J^ AGeV [21]. The small 
asymmetric contribution in the first system is neglected in evaluating mean values 
and variances. 

slightly larger than the average over all fission fragments produced. We as- 
sume that*** fission proceeds from a mean parent nucleus of mass A = 220 
with z/ = 6 the averaged number of emitted neutrons from the symmetric and 
asymmetric modes. On tabic 2 the cross sections, mean values of A and Z, 
variances and local variances arc presented and compared to the results of 
symmetric fission in (U + Pb) [21]. 



4.-5 Kinetic energies 
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Fig. 17. Mean kinetic energy of fragments from 1 A GeV ^^^U + p fission as a 
function of their atomic number Z. 

The average kinetic energy of fission fragments as a function of their atomic 
number is given by Ekin (Z) = i.mo A{Z)Vf where A{Z) is the average mass 
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Table 3 

The mean kinetic energy values of single fission fragments as a function of the 
atomic number Z. 

number obtained from the measured isotopic distribution, and V/ is the cor- 
responding fragment velocity as shown in Fig. 9a. The mean kinetic energies 
are given in table 3 and shown in Fig. 17. They decrease smoothly for higher 
atomic numbers and show a broad maximum around Z = 32. The mean value 
of the kinetic energy, calculated by using the Z-distribution of Fig. 14a gives 
(76 ± 3) MeV, corresponding to a mean total kinetic energy release of (152 ± 
6) MeV. 

J!i..6 Comparison with simulation codes 

The set of isotopic distributions of fission fragments provides a real challenge 
for simulation codes. A fast and simple formulation developed by Silberberg 
and Tsao [51] is shown in Fig. 18 and 19. The model reproduces the orders of 
magnitude of the isotopic distributions, but the asymmetric fission component 
is overestimated in magnitude and symmetric components are too narrow. 
Isotopes of elements at and below Z = 35 are poorly reproduced. A jump at 
Z = 35 in the code produces the unrealistic gap in the mass distribution at A 
= 90. 

The second code widely used is the LAHET-code (Los Alamos [24] or Oak 
Ridge version) [52]. The agreement with this code is not satisfying. In Fig. 
18 and 19 the Oak Ridge version has been used for comparison. For neutron- 
deficient isotopes, cross sections are overestimated, and symmetric fission is 
underestimated. An odd-even staggering is produced not seen in the experi- 
mental data. 

Parallel to our experimental program different approaches to new simulation 
codes are published or are under way. The first reaction step in spallation. 
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Fig. 18. Comparison of measured isotopic distributions (full points) to calculations 
with the formula of Silberberg et al. [51] (crosses), by LAHET [52] (triangles), and 
by our improved codes [25,26,27,28] (diamonds) for a selection of elements 

the intranuclear cascade induced by the primary proton collision, is simulated 
by the new version of the code INCL [28]. At this stage, a distribution of 
prefragment nuclei are predicted with an associated excitation energy and 
intrinsic spin. This is then the physical input to a system of codes simulating 
the second step of spallation, the de-excitation phase. New results on the 
physics of particle evaporation were introduced in the statistical de-excitation 
code ABLA [25,26] and on fission in the fission code PROFI [27]. 

The fission code PROFI is a semi-empirical Monte-Carlo code which calculates 
the nuclide distribution of fission fragments. It is theoretically based on the 
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Fig. 19. Mass distribution of fission fragments as a function of mass number com- 
pared to simulations; Silberberg [51] (crosses), LAHET [52] (triangle) and to our 
improved codes [25,26,27,28] (diamonds). Experimental data are shown as full sym- 
bols with error-bars. 

application of the statistical model of nuclear reactions to the concept of fission 
channels. In this model, the population of the fission channels is assumed to 
be basically determined by the number of available transition states above the 
potential energy surface near the fission barrier. Several properties, however, 
are finally determined at the scission time. A full description of the model is 
given in [27]. 

The barrier as a function of mass asymmetry is defined by three components. 
The first is the symmetric component defined by the liquid -drop potential by 
means of a parabolic function with a curvature obtained from experimental 
data [53]. This parabola is assumed to be modulated by two neutron shells, 
located at mass asymmetries corresponding to neutron numbers N = 82 (spher- 
ical neutron shell) and N = 90 (deformed neutron shell). We assume that the 
mass-asymmetry degree of freedom at the fission barrier is on the average 
uniquely related to the neutron number of the fragments. The shells are rep- 
resented by Gaussian functions. These shells are associated with the fission 
channels Standard I and Standard II, respectively [54], while the liquid-drop 
potential is associated with the symmetric fission channel. The population of 
the fission channels is proportional to the level density around the correspond- 
ing dips in the potential at saddle at the available excitation energy. Shells 
are washed out with excitation energy [55]. The heights and the widths of the 
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Gaussian curves representing the shell effects and additional fluctuations in 
mass asymmetry acquired from saddle to scission, are derived from experimen- 
tal data [27]. The mean values of the neutron-to-proton ratio for the channels 
Standard 1 and Standard II are deduced from measured nuclide distributions 
after electromagnetic induced fission of ^ss^j po]. The charge polarisation for 
the symmetric fission channel and the fiuctuations in the neutron-to-proton 
ratio for all channels are also considered by describing the potential in this 
degree of freedom again by a parabolic function [48] . Assuming that the equi- 
libration in this variable is fast compared to the saddle-to-scission time, this 
potential was calculated in the scission configuration. Since the shell effects of 
the nascent fragments at scission, which are strongly deformed on the average, 
are not known experimentally, only macroscopic properties are included in the 
calculation of the charge polarisation of the symmetric fission channel. 

Consequently, the two fission pre-fragments are obtained. Their excitation en- 
ergies are calculated from the excitation and deformation energy of the fission- 
ing system at the scission point. The fission probability and the consecutive 
de-excitation of the fission fragments via particle evaporation are conducted 
by the routines of the ABLA code [25,26]. Coupling to the cascade code INCL 
[28] gives a reasonable reproduction of our data, as demonstrated in Fig. 18 
and 19 (diamonds). Except for an excess of asymmetric fission the simulation 
meets the experiment. Either INCL provides a too small excitation energy 
for the most peripheral collisions, or the excitation energy in the statistical 
de-excitation chain is underestimated or, finally the parameter set in PROFI 
has to be reajusted in order to reduce the share of asymmetric fission. 

Jf..! Production of very neutron rich isotopes 

The experimental method of 1 A GeV U-projectile fission opens a very effi- 
cient way to produce secondary beams of neutron-rich isotopes. In order to 
provide experimental cross sections for such a project, the production of very 
neutron-rich fragments by 1 A GeV U-|-p collisions is evaluated renormalising 
the values obtained from a previous dedicated measurement for U-|-Be [8] . The 
production of fission fragments differs in both systems because of a relatively 
larger contribution of asymmetric fission with U -|- Be. Normalising the el- 
emental cross sections for U + Be to the corresponding element for U + p, 
the production of very neutron-rich can be deduced reliably down to a level 
of 1 nb. Not only the falls of the isotopic production with increasing neutron 
excess are the same for U -|- Pb and U -|- Be [20,8], but also the production 
is found to be the same in the gap of overlapping masses presently covered in 
U + p. 

Fig. 20 illustrates four examples of isotopic distributions. Such exotic beams 
in spite of their low intensities will be of importance to test predictions for 
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Fig. 20. Comparison of cross-sections from ^ss^j ^ ^ GeV + p (points) to ^^^U 0.75 
AGeV + Be [8] (triangles). These last data are renormalised in a{Z) to the present 
elemental cross sections. 

the vanishing of shell-effects at N = 82 and 50, and the possible appearance 
of neutron halos in the region of tin and nickel isotopes with large neutron 
excess. 



4-8 Parent fissioning nuclei 

The reconstitution of the charge Zq and mass Aq of the parent fissioning 
nuclei is of importance to constrain critical parameters of the intranuclear 
cascade phase: the transfer of excitation energy and of momentum to the pre- 
fragments, the value of the density of intermediate states and of the viscosity 
of highly excited heavy nuclei. The present experimental findings should con- 
verge with all other results of ^ss-q x a GeV + p [17,44,46] on an unified 
description of a first stage of nucleon-nucleus collision described by a cascade 
followed by an evaporation stage among which fission occurs. 

What is learnt about the fissioning parent nuclei in the present work is listed 
below: 

• Isotopic distributions lead to a mean fission parent nucleus ^^°Th. 
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• For a given element the increase of the mean velocities with the mass of the 
isotopes (Fig. 8) indicates that parent elements in the range of Zq < 90 do 
contribute to the fission. 

• Velocities of fragments Fig. 9a show that the region of 50< Z < 64 arises 
from parent elements of ~ 90, while the region of 30 < Z < 37 indicates 
lower fissioning elements at Zq ~ 84. 

• The variance of the fission fragment velocities ayf is extracted from the 
external slopes of the observed velocity distributions and presented in Fig. 
9b. The value of 0.13 cm/ns, the same for all fission-fragment elements, is 
correlated to the variance of the recoil momenta generated in the nuclear 
cascade for the ensemble of the fissioning parent nuclei. At AA = 18, the 
analysis of spallation presented by J. Taieb [17] gives ap — 350 MeV/c as 
the variance of the recoil momenta of the central parent nucleus, ^^°Th. This 
variance is transferred to the fission fragments and leads to a contribution 
of 0.09 cm/ns in the fissioning system. The ensemble of parent elements 
and isotopes contributes also by the fiuctuations of their TKE-values. A 
pair of fission-fragments issued from a fixed parent nucleus shows a vari- 
ance of its velocities. Taking from thermal-neutron induced fission of 235]j 
the measured value of (Ttke = 6.0 MeV as a lower estimate, a contribution 
of 0.02 cm/ns to the variance of fission fragment velocities is calculated. A 
much larger contribution is generated by the difference in TKE-values of 
the many nuclei present. The variance in the Zo-distribution of the fission- 
ing parent nuclei of a^g/Zo of 4.6% is taken from the measurement of B. 
Jurado [46]. Including the dispersion due to a spread over several isotopes 
in each element, and converting TKE-values into fragment-velocities, the 
contribution of the fissioning ensemble is evaluated to 0.09 cm/ns. The lo- 
cation straggling in the thin H2-target induces a negligible contribution of 
0.01 cm/ns to this variance. Adding quadratically all contributions, a total 
variance of 0.13 cm/ns is obtained which coincides with the value extracted 
from the observed velocity distributions. 

• The mean recoil velocities of the parent fissioning nuclei arise from the first 
phase of the reaction and contribute to the momentum losses measured [17]. 
The transferred momentum of -150 MeV/c measured for the evaporation 
residues at A A = 18 should be close to the momentum of the mean fissioning 
parent ^^°Th. This recoil momentum is equivalent to a recoil velocity of - 
0.04 cm/ns, indicated by a line in Fig. 9c. The recoil velocities presented in 
Fig. 9c, in the range of -(0.04-0.13) cm/ns are on the average larger by a 
factor of 2 and show a trend towards larger values for fighter elements. 



4-9 Comparison of system U + p and Pb + p 

Fig. 21 and Table 4 demonstrate the basic difference between the collision 
systems Pb -|- p and U -|- p. Both reactions have nearly the same reaction 
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Fig. 21. Comparison of cross-sections for evaporation and fission residues as a 
function of the mass-loss AA for the two coUision systems Pb -|- p (triangles) and 
\J + p (full points) at 1 A GeV. 



reaction 


(barn) 


CTfis 

(barn) 


O'EVR 
(barn) 


U + p 
Pb + p 


1.99± 0.17 
1.84± 0.23 


1.53 ± 0.13 
0.16 ± 0.07 


0.46±0.08 
1.68±0.22 



Table 4 

Comparison of reaction cross sections of 1 A GeV ^^^U -|- p to 1 A GeV 
208pb + p [15]. 

cross section, 1.84 b for Pb + p and 1.99 b for U + p. Fragmentation is the 
main reaction channel for Pb+p, with a = (1.68 ± 0.2) b, whereas fission 
dominates the U-|-p reaction with a cross section a = (1.53 ± 0.15) b. 

Averaging over all cross sections from evaporation and fission residues, alto- 
gether a mean neutron number of 127 ± 1 is found for reaction products of 
238u -|- p at 1 A GeV, i.e. 19 ± 1 neutrons are liberated in the reaction. This 
number does not include neutrons bound in the lightest elements ( Z < 7 ), 
and should be an upper limit for the number of neutrons produced in thin 
targets. It compares with the neutron multiplicity of 20 measured with 1.2 
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GeV proton on depleted uranium [56] when extrapolated to thin targets. 



5 Conclusions 



We have measured isotopic production cross-sections for about 733 fission 
fragments produced in 1 A GeV ^ss^j _j_ p collisions with a mean accuracy of 
10 %, down to values of 20 fib for neutron-rich and 0.1 mb for neutron-deficient 
isotopes. The population of isotopes resulting from fission is pictured in Fig. 
22 on the chart of nuclei. The isobaric slope of the production cross section 
towards the neutron rich side of the chart of isotopes is found to be the same 
for hydrogen as for heavier targets investigated earlier. 




Fig. 22. Two-dimensional plot of the isotopic cross sections for fission-fragments 
obtained in 1 A GeV ^^^U -|- p shown on the chart of isotopes with squares indicating 
stable isotopes. Colours correspond to increasing cross sections according to the 
logarithmic scale indicated. 

The fission cross section, including the production of very light elements, 
amounts to (1.53 ± 0.15) b which compares well with previous results [35,45] 
as with the recent result of B. Jurado [46]. The total cross section of fission 
plus evaporation residues amounts to (1.99 ± 0.17) b, close to the calculated 
geometrical cross section of 1.96 b [42,47] and to the INCL calculation 1.94 b 
[28]. 
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Our results show a good quantitative agreement with previous isotopic dis- 
tributions of cross sections measured for rubidium and cesium with on-hne 
mass-separator techniques except for a discrepancy in cross sections for the 
6 most neutron-rich isotopes of cesium [3]. The coherence of the present ex- 
perimental method validates our mass-calibration relatively to all other iso- 
topes. Parameters characterising fission and its variances are extracted. The 
elemental distribution and the kinetic energies of the reaction products are 
of importance concerning chemical corrosion and the composition of nuclear 
waste produced by spallation of actinides. 

The analysis of the momentum and velocity spectra of fission- fragments as 
function of A and Z gives an indication on the nature and the domain of 
excitation energies of the fragments arising from the first step of the (U -|- p)- 
interaction and undergoing fission. Comparison with previous measurements 
shows that the symmetric fission cross section (1.42±0.20) b is almost the same 
for U -I- p as for U -|- Pb, (1.4 ± 0.2) b. In this last case the first phase of violent 
abrasion leads to a wider range of more excited and lighter fragments. The 
fission probabilities of these fragments are smaller than for actinides close to 
the U-projectiles which are mainly produced in (U -|- p)-collisions. Finally the 
cross section to observe symmetric fission remains the same for both systems. 

There is certainly not a perfect overlap between the distribution of fragments 
identified as evaporation-residues and the fissioning parent nuclei. The par- 
ent nuclei close to U have high fissilities. For elements of decreasing Z num- 
bers, fissile isotopes become more and more neutron-deficient. Compared to 
the evaporation-residues the fissioning-parent nuclei are shifted towards the 
neutron-deficient side. They are born from prefragments of higher excitation 
energies and their recoil momenta may be higher than for those ending as 
evaporation residues. 

Measured isotopic and mass distributions are compared to results obtained 
with simulation codes commonly used. Large discrepancies are shown, and the 
need for improved reaction models is underlined. Models developed recently in 
the collaboration are discussed [25,26,27,28] and the improvements achieved 
are demonstrated. 
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7 Appendix 1 

1.1 Demonstration of equations (7) and (8) 

The transformation from the c.m.-system to the lab-system is given by: 



/^f sin 09 , , 

7tan0^ 7 11 

with ± signs for forward and backward angles ip^ and respectively. To 
transform back to the c.m.-system for the limiting case (pum = a, we solve 
equ. (11) for this case: 



7tana=— — 12 
K ± cos ip 

with K = /9o / I3f- It is convenient to solve equ. (12) as a quadratic equation 
in tan(<^/2), which has 2 solutions corresponding to tp^ and (p^ given for 0;^^ 
= q; in the lab-system. 



tan((^^'^/2) = tan 8/2 (13) 

with an auxiliary variable 



sin 5 = 7V K"^ — 1 tan a (14) 

Equ. 13 shows that the angles (p^ and</?^ depend via eq. (14) on all 3 variables 
a, Po and Pf. The transmission of each of the cones may be evaluated with the 
total transmission obtained by summing T^ and Tq. This method we applied 
in our evaluations until now [13,15,16]. However, introducing the two solutions 
of equ. (13) into equ. (3) we obtain another very compact expression for Tf^: 



cos 5 

Ta^l- — ^- ^ 15 

1 + (7tanQ;)^ 



Equ. (15) depends via the auxiliary variable 6 on the fission velocity. To il- 
luminate the variable 5 we introduce (ps — {(p^ + <p^)/'2, and (pA — {^p^ - 
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(p^)/2, or (p^ = (pY, + <^A and (p^ = (p^,- ^a- Equ. (13) transforms to the new 
variables using the trigonometrical sum relations for tan (/7s and tan(/?A- 



tan ifs — I = tan 6 tan (/?a — 7 tan a (16) 



Angle (pA depends only on two variables, 7 and a, which represent the rela- 
tivistic beam and the FRS spectrometer. The fission velocity /3/ enters into 
the angle (p^. Between (p^'^ and </7s,a another trigonometrical relation holds 

cos (p^ + cos (p^ —2 cos (fs cos (fA 

This relation introduced in equ. (3) gives for the transmission : 

7n = 1 - cos (fs cos (pA (17) 
Comparing equ. (15) and (18) a relation between (/?e,a and 5 follows: 



coso = = cos(p-£\ 1 + (7tanQ;j^ 

COS</?A 

The three expressions equ. (3), (15) and (18) are exact, identical representa- 
tions in different angular coordinates. All expressions derived for the transmis- 
sion until now depend on f3f. The conversion of (3f into the measured quantity 
use of equ. (6) and (14) : 



sin 5 = ^ ° S tan a = (2 - T^) (18) 



Also equ. (6) and (19) are exact without any approximations. 

For our experiment with 7 = 2.07 and a — a — 14.8 mrad, is larger 
than 5 by only 0.1 mrad. It is a good approximation to replace </7e, the mean 
value of (p^ and (p^, by the auxiliary variable 5. A next approximation for the 
transmission follows: 



7b = l-cos(/?E (19) 

which deviates from the exact value by less than 0.2 % in the worst case of 
Z = 28. Neglecting terms in and useful approximations derived from 
eqs (16) and (17) are: 
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tpY, = 8 and <^A = 7« 

— ^i: -\- 70; and — ^y, - loc 
Prom equ. (20) and (6) we derive equ. (7) and (8). 
With I3iijn = a — 1 follows: 



aV7^ - 1 Plim .„„s 

tan(/?s/2 = — -5 = ^^3— (20) 

Finally with equ. (21) Pf and Tq can be presented as a function of (3app, see 
equ. (7) and (8) given in section 2.5. 



1 + mpp/(3um? 

and 

Pf — Papp[i + iPlim/'^Papp) ] 

q. e. d. 
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Table 5: Fission fragment isotopic cross sections mea- 
sured in the present work. Statistical uncertainties are 
given on the last significant numbers, excluding the 10% 
systematical uncertainties. 
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